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Abstract 

We study the merging of two independent Bose-Einstein condensates with arbitrary initial phase 
difference, in the framework of a one dimensional time-dependent Gross-Pitaevskii model. The 
role of the initial phase difference in the process is discussed, and various types of phase-sensitive 
excitations are identified. 
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I. INTRODUCTION 



A large class of applications of Bose-Einstein condensation of dilute gases involves one way 
or another, the merging of two initially separate condensates. Matter wave interferometry 
[l| and the quest for a continuous atom laser j^jlsl, are only two prominent such applications. 
The process of merging requires the controlled unification of two trapped clouds. Prepared 
independently in separate traps, the two clouds are expected to be united into a single, 
trapped condensate of a well defined state. One example of the successful merging of two 
independent condensates has been demonstrated by Chikkatur et al. js], through the use of 
optical tweezers to transport one condensate from their production chamber to merge it with 
a previously prepared one. As expected, the final cloud was found to contain more atoms 
than either of the two initial clouds, but less than their sum; presumably due to evaporative 
losses. 

Theoretical investigations and modeling of the process have so far been rather limited, 
with a number of relevant questions still open. More precisely, in related work 0-3, the 
influence of an initial phase difference between the two condensates on the dynamics of the 
merging, the question of the phase of the final condensate, as well as the type of phase- 
sensitive excitations were not addressed. These aspects are of vital interest as attested 
by their frequent recurrence in the relevant experimental literature . Given that two 
independently formed condensates will have a random phase difference, it would seem that 
issues pertaining to the phase merit attention not only from a fundamental point of view, but 
also for the applications envisioned. A case in point is highlighted by the recent experiments 
of Jo et al. [TJ reporting the dependence of heating and atom loss during the merging on 
the phase difference of the fragments. Motivated by the above theoretical and experimental 
developments, our aim in this paper is to explore somewhat further these aspects. Although 
the work in this paper is limited to zero temperature, it does provide useful insight on the 
role of the initial relative phase on the merging. 

II. THE SYSTEM 

The system under consideration pertains to two independent elongated condensates (L 
and R) consisting of a large number of bosonic atoms cooled into the lowest eigenmode of the 
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corresponding harmonic trap. The merging of the two condensates is achieved by bringing 
the two traps together in a controlled and adiabatic manner. 



A. Model of the merging 

Following js], we have investigated this process in the context of a one-dimensional model. 
As we will see later on, despite its simplicity this model is capable of capturing many of 
the phenomena that take place during the merging process, and have not been addressed in 
earlier related theoretical work js, 6|. 

To be consistent with , as well as the experimental setup for condensate merging 
[sl, we assume two nearly identical harmonic traps with confining frequency 00. As the two 
traps move towards each other, the global potential experienced by the trapped atoms can 
be modeled by a double- well potential of the form js, 9| 



V{x,t) = ^i\x\-s{t).lf 



(1) 



in dimensionless units. The function s{t) determines the details of the merging (i.e., speed 
and time scale Tm), which has to be adiabatic so that any kind of excitations due to movement 
of the traps are suppressed. To this end, the transport of the condensates must take place on 
a time scale much larger than the characteristic time-scale of excitations along the merging 
direction, as well as the time-scale of interatomic interactions js, 6|. Moreover, excitations 
can be minimized by appropriately choosing the profile of s{t). 
Throughout our simulations, s{t) has been chosen as 



s{t) 



cos 



1 for t = 0, 

for < t < Tr, 
for t > T^. 



(2) 



In the beginning of the merging (i.e., at t = 0) we have two well-separated traps, and 
the potential ([1]) exhibits minima at x = ±/, with / chosen sufficiently large. During the 
merging, i.e., for < t < the two traps approach each other, and the two minima of the 
double-well potential are located at x = ±s{t)l. Accordingly, the barrier between the two 
wells also decreases in this regime and in the end of the merging (i.e., at t = Tm) we have 
complete overlap of the two condensates. After the merging, i.e., for t > T^, the potential 
remains a single harmonic well. 



We choose to describe the evolution of the system by the time-dependent Gross-Pitaevskii 
equation (GPE), which in dimensionless units reads 

d 1 

(3) 

with the function normahzed to unity. The nonhnearity parameter g is proportional 

to the total number of atoms and the corresponding s-wave atomic scattering length. We 
use harmonic oscillator (h.o.) units, i.e. -y/^, and hhj for length, time and energy 
respectively, where m is the atomic mass. The energy of a solution of Eq. ([3]) is given by 



= y"dx Q|V^(x,t)p + \/(x,t)|^(x,t)p 



= -Ekin + Epot + -Eint, (4) 

and is a conserved quantity provided that the external potential does not depend on time 



10|. The lowest energy solutions of Eq. 02]) to a given potential V{x) can be written as 

= ^(a;)e-'^*/'"'+''^. (5) 

Here, if){x) is a real function whose squared modulus is the atomic density, while is a 
random phase which emerges, as a broken gauge symmetry, in the process of creating a 
condensate 10|. The chemical potential /i is given by 

/i= (^km + ^pot + 2^int). (6) 

In general, the time-evolution of a condensate which is initially in its ground state, is uniquely 
determined by the two parameters fi and (f. 

As mentioned earlier, in the beginning of the merging, the two identical condensates 
are considered to be independent. This means that each condensate experiences its local 
harmonic potential only and thus, without loss of generality, the initial state of the system 
is given by 

t = 0) = <I>l(x) + <l>R(x)e'^'^^", (7) 

with $L(R)(3;) the lowest-energy solution for the isolated left(right) harmonic potential re- 
spectively, and Ay9in = fR — fh the initial phase difference between the two condensates, 
where (/3r(l) is the phase of the right (left) condensate. 



FIG. 1: Initial conditions used in the numerical simulations. Symbols: atomic densities |^'(x,0)p 
corresponding to the ground state of Eq. Q for the double- well potential V{x,0). Lines: atomic 
density distribution for two identical independent condensates trapped in harmonic potentials 
centered at a; = ±4 (see Eq. [7]). Harmonic oscillator units are used for x. 

In principle, however, this is not the case. According to Eqs. ([T]l2]) the two condensates 
at t = experience the double-well potential V{x,0) = {\x\ — /)^/2. Hence, to ensure inde- 
pendence of the two condensates, throughout our simulations we had to choose sufficiently 
large separation /. In that case, the lowest-energy solution of the system \l/(x,0), is well 
approximated by Eq. (see Fig. [1]), but the two condensates are always in phase, i.e., 
Ay^in = . We had therefore to introduce by hand the initial phase difference between the 
two condensates, by multiplying ^'(x, 0) by a factor e''^''''", for x > 0. Although the resulting 
state is no longer longer a ground state of the double-well potential, one may readily check 
that the increase in energy is very small due to the small value of |\I^(x, 0)p at x = 0. Finally, 
it is worth noting that the symmetry of the model for alH > ensures that the initial phase 
difference, is not affected by the movement of the traps. 



III. NUMERICAL SIMULATIONS 



We solve Eq. ([3]) using a time-splitting spectral method on a spatial grid of 2048 
points ranging over 20 harmonic oscillator units. The size of the applied time step was 
10"'^ and the numerical performance - estimated by the conservation of particle number and 
energy - was found to be very good. Simulations were performed for various initial phase 
differences and merging times, in the absence, as well as in the presence, of interatomic 
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interactions. In the latter case, the dimensionless nonhnearity parameter was chosen to be 
g ^ 8.5, which characterizes an intermediate regime of interaction strength. Throughout 
our simulations we have investigated the dynamics of the merging from the perspective of 
two quantities, namely the energy of the system and the phase of the final condensate. 

It is worth keeping in mind for the following discussion that the system under investigation 
is invariant to changes of the global phase. This is obvious from the fact that the GPE 
which governs the evolution of the total wavefunction of the system remains invariant under 
transformations of the form t) — )■ e*-^\Ef (x, t), where x is a constant global phase. Hence, 
such changes do not affect quantities that are determined only by densities |\E'(a;,t)p and 
phase differences (such as energy). Moreover, this invariance gives us some freedom in 
choosing a reference phase. Throughout our simulations, the reference phase has been 
chosen as (/:>l = [see Eq. ([7j)]. 

A. Energy 

The time evolution of the energy of the system, as determined by Eqs. (jSH]), for various 
initial phase differences is depicted in Figs. [2](a) and (b). At t = we find -Eid(O) = 0.5 
and -Eint(O) = 1.20 for ideal {g = 0) and interacting gas {g ~ 8.5), respectively. Clearly, due 
to the large initial separation of the traps, there is no remarkable dependence on the initial 
phase difference. At the very early stage of the merging i.e., for < t < 10, the overlap 
between the two condensates is negligible, and thus the energy remains very close to its 
initial value; a fact which also confirms the initial independence of the two condensates as 
well as the adiabatic nature of the merging process for the chosen parameters. For longer 
times (i.e., for 10 < t < T^), the condensates' wavefunctions begin overlapping in space, 
and the energy of the system depends crucially on the initial phase difference. In particular, 
as a general observation we note that the energy of the system increases as we increase the 
initial phase difference A(/?in. At t ~ 15, we observe a local minimum in energy for small 
values of A(p[^. It occurs for both ideal and interacting gases; albeit at slightly different 
time instants. Hence, its appearance is not associated with the nonlinear character of the 
GPE, but is a feature of the particular double-well potential ([T]) and the energies of its two 
lowest eigenstates j^. 

In an ideal scenario, where all types of excitations are suppressed, one expects the final 
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FIG. 2: (Color online) (a) Time-evolution of the total energy of an interacting gas during a merging 
process, for various initial phase differences, (b) The same for an ideal gas. (c) Excess energy at 
the end of the merging as a function of the initial phase difference Ay^in. The dashed line through 
the crosses is the theoretical result for ideal gas, given by Eq. Q. For symmetry reasons, all 
the plots can be extended to A(/?in < by the transformation Ac^in — > — A(^in. Merging times: 
Tjn = 22.5 (a,b) and Tm = 23 (c). Energy and time are in units of huj and f^^"*^, respectively. 
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condensate to be in the ground-state of the harmonic potential V{x,T^), with the corre- 
sponding energy being -Ej^l^'* = 0.5 and -Ej^jft^ = 1.76, for ideal and interacting gas respec- 
tively. In our simulations, however, although we have ensured adiabaticity, the final energy 
Eid{int)(Tni) of an ideal(interacting) gas may exceed the corresponding ground-state energy 

^idjnt)- Fig- Plo* excess energy in the system t/id(int) = ^id{int)(7'm) -^JlJ^t)' 

a function of the initial phase difference. In the presence of interactions, the excess energy is 
negligible for Aipi^ < 7r/4 and, to good accuracy, the system is in its ground state at the end 
of the adiabatic merging. For an ideal gas, the corresponding regime is much narrower as it 
pertains to very small phase differences Ay^in < ir/S. Moreover, in both cases we observe a 
rapid increase of the excess energy for increasing A(fi^, but for different reasons. 

In the case of a noninteracting gas, the problem under consideration reduces to the 
problem of a single particle in a time-varying double-well potential, because Eq. ([3]) reduces 
to the Schrodinger equation. Under the assumption of adiabatic merging, one may readily 
obtain analytic expressions for the occupation probabilities of the ground state (po) and the 
first excited state (pi) of the final harmonic potential V{x,Tya)- The population imbalance 



is given by [12 1 



PO - J>1 = C0s(Av9in), (8) 

and thus the excess energy (in units of u) reads 

t/id = [1 - cos(A^i,)] /2. (9) 

This theoretical curve is also drawn in Fig. M^c) and shows very good agreement with our 
numerical results. Thus, in the case of an ideal gas, at the end of the merging the system is 
in a superposition of the ground state and the first excited state of the harmonic potential 
V{x,Tj^). In the extreme case of A(pi^ = n, the final state is basically the first excited state 
of the harmonic oscillator, i.e., the first asymmetric Hermite polynomial, which has a node 
at X = 0. Hence, the final density exhibits two distinct peaks, while it vanishes at the center 
of the trap. 

In the case of an interacting gas, the excess energy observed in our simulations is due 
to a dark soliton which is formed adiabatically during the merging process. The initial 
phase difference determines the depth of the soliton, and thus its dynamics. As before, 
merging is impossible for Aifi^ = vr, where we have the formation of a static black soliton 
and the final density in the single well exhibits two distinct peaks [see Fig. [3](f)]. For smaller 
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phase differences, however, we have the formation of shallower solitons [Figs. [3]^c-e)], which 
oscillate back and forth in the trap [see Fig. |3l^a)]. This is in agreement with known results 
of soliton dynamics 13|, and similar to the formation of vortices in the merging of three 
Bose-Einstein condensates with appropriate phase differences fl^ . 

We have seen therefore that, for both ideal and interacting gases, the adiabatic merging 
results in a condensate whose density exhibits a dip. The origin of the dip, however, is 
fundamentally different in the two cases, and this fact is expected to be reflected in the 
dependence of the dip's depth on the initial phase difference. In Fig. [3]^b) we plot the the 
dip's depth at x = (i.e., n = |\l/(0,t)p for t > Tm), normalized to its value Uq estimated 
for A<y5in = 0, as a function of A(pi^. Using Eq. ([H]), in the absence of interactions we obtain 
n/rio = [l + cos(A(/9in)]/2. For an interacting gas, however, the ratio n/uQ behaves differently 
for varying A(/9in. In particular, we see that the dip remains rather shallow for A(^in < 7r/2 
and becomes deeper rather abruptly as we increase Aipi^i further. In any case, our results 
show that one may get a rough estimate of the relative phase between two condensates by 
looking at the density profile after their adiabatic merging. 

In closing, we would like to point out that Figs. [2]and[3]^b) can be extended to negative 
phase differences, taking the mirror images of the corresponding curves for Ay^in > 0, with 
respect to the vertical axis. This is due to the symmetry of the system under exchange of 
the two wells, which together with the global-phase invariance, imply that quantities that 
depend only on the density and the phase difference, are symmetric under the transformation 
Ay^in — 7- — A(y9in. Indeed, the system is initially prepared in the state with real $L{R)(a;)- 
As discussed earlier, the dynamics of the system, from the point of view of densities and 
energies, do not change if we multiply the entire state by e~^^'^''\ Doing so, the initial 
condition reads t = 0) = $L(2;)e~'^'^'" + $r(x). Given that the symmetry of the 

system with respect to the two wells is preserved throughout the merging process, we can 
exchange the two initial wells obtaining t = 0) = ^h{x) + $Pi,(x)e~'^'^'", which differs 
from the initial condition by the sign of the phase difference. 
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g = 8.5 □ 
■^(l + cos(A^i„)) 








FIG. 3: (a) Time-evolution of the atomic density distribution in arbitrary units after 

the merging of two identical condensates with At^in = 37r/4. Harmonic oscillator units are used 
for space and time, (b) Relative depth of the dip in the final atomic density distribution at the 
center of the trap, for various values of At^in- (c-f) The entire distribution is depicted only for 
Aifia = 0, 7r/3, 27r/3, and tt. 

B. Center of Mass Motion 



To gain further insight into the nature of the excitations that occur during the merging, 
it is interesting to investigate the dynamics of the center of mass motion 10 1 



x{t) 



dx\'^{x,t)\'^x. 



(10) 
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The evolution of x{t) as a function of time for various initial phase differences, is depicted 
in Fig. m For all the choices of parameters, the behavior of x{t) was found to be well 
approximated by 

x{t) = Csm[{t-Tj + e], (11) 

with the amplitude C and the phase 6 depending on the initial phase difference and the 
strength of interatomic interactions only. 

As depicted in Fig. IH^c), for g = the amplitude of the oscillations varies sinusoidally 
with Aiy^in, and is well approximated by 

C ~ 0.75 sin(A(/?in). 

The presence of interactions deforms this symmetric behavior around 7r/2, since C attains 
its maximum for higher values of Aipij^. As far as the phase 6 is concerned, for an ideal gas it 
depends on only, while for (7 > it also acquires a dependence on Aipi^. According to our 
simulations, the excitations discussed in Sec. IIII Al (i.e.. the presence of a soliton for g ^ 0, 
or the population of the first excited state for g = 0), cannot describe all aspects of the 
center of mass motion. The above observations, as well as additional details not described 
here, show that besides the excitations discussed earlier, the final condensate performs dipole 
oscillations with an amplitude depending on the initial phase difference. 



C. Final Phase 



or an easy modification thereof. 
, Q], however, suggest a fast decay of phase-sensitive 



As we saw in the previous subsections, for an interacting gas at the end of an adiabatic 
merging, we can distinguish between two types of phase- sensitive excitations in the final 
condensate, namely dipole and soliton-like excitations. In the one-dimensional model under 
consideration, such excitations do not decay, and thus the state of the final condensate 
cannot be expressed in the form of ([5j 

Recent experimental observations 
excitations in a three-dimensional merging setup, which results in an increase of the tem- 
perature on the order of hw/k^ ~ nK, where ke is the Boltzmann constant. This is in 
agreement with our estimates for the excess energy in the final condensate [see Fig. [21(c)]. 
This amount of energy can be dissipated by evaporative cooling, with a small loss in the 
number of atoms [sl,!^]. Most importantly, such a cooling mechanism is not expected to 

11 




FIG. 4: (Color online) (a) Displacement of the center of mass in h.o. units with no interatomic 
interactions for Tm = 25lo~^. (b) Same as (a) but with interactions g = 8.5. (c) Amplitude of the 
dipole oscillation in dependence of the initial phase difference Aipi^. 

destroy the phase information carried by the final condensate; a very crucial issue for many 
applications (e.g., interferometry atom lasers {2, 3|). 

As an attempt to trace the final result of such a dissipation mechanism after the merging 
(i.e., for t > Tm), we propagated Eq. ([3]) in imaginary time (i.e., replacing idt by dt). This 
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FIG. 5: (a) The phase of the final condensate (pi, as a function of the initial phase difference 
Aifii-i in the noninteracting case. Merging times are Tm = 21.52lo~^ (crosses) and Tm = 26.52uj~^ 
(squares), (b) Same as (a) but with interactions and merging times Tm = 23u)^^ (crosses) and 
Tin = 23 Auj~^ (squares). The line is Eq. (fT2]) applied to the upper curve. 

imaginary-time evolution, does not affect the complex phase of \l/(x,t), but only the density 
|\l/(x, t)p and appears as particle loss. Thus the function \I^(x,t) has to be continuously 
renormalized and ends up eventually in a stationary state of the form ([5]), enabling us to 
extract the chemical potential and the phase ipi of the final condensate. 

According to our simulations, the final phase (pf depends mainly on the merging time 
Tm, the initial phase difference Ay^in, and the interaction strength g. More interestingly, 
our simulations reveal the distinct roles of these parameters on the the phase of the final 
condensate. In Fig. [5l we present results for ideal and interacting gases and for various 
merging times. These results pertain to the initial state ([7]) where, without loss of generality, 
we have chosen as a reference phase (pi^ = and thus effectively A(fi^ = ipn. 

In the absence of interactions, and for a given merging time, the final phase depends 
linearly on the initial phase difference Ay^in, with the corresponding slope being approxi- 
mately equal to 0.5 [see Fig. |5]^a)]. For increasing merging times, the entire curve shifts 
upwards linearly with Tm, without any noticeable effect on its form. Hence we have 
ipf ~ /3(0,Tin,0) + A(/)in/2, where f3{g,T^,ipi^) is a function of g, Tm, and ipi^. In view 
of the global-phase invariance discussed in Sec. Illlt the same results hold for arbitrary val- 
ues of the reference phase ipi^ with f3{g, T^, ipi^) = f3{g, T^, 0) + ipi^. Hence, the phase of the 
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FIG. 6: The function a in dependence of the merging time for interacting (squares) and non 
interacting (crosses) gases. 



final condensate is well approximated by 



(/?£ = « + 



(12) 



where for the sake of simplicity we have set ai^g^T^ = /3((7, Tm, 0), with a being a linear 
function of the merging time (see Fig. |6]) 16|. Note that Eq. ( IT2l) is invariant under the 
exchange of the two wells in agreement with the symmetries discussed in Sec. Illll 

For an interacting gas (i.e., for g ^ 0),we have a similar behavior of the final phase as we 
vary all the relevant parameters. As depicted in Fig. ^h), for |A(/?in| close to zero the final 
phase is well approximated by Eq. ( IT2|) . The parameter a^g^T^^) is a linearly increasing 
function of the merging time, with the slope being determined by the interaction strength g 
(see Fig. [6]). For A(^in ^ ±7r, however, we observe strong deviations of the final phase from 
Eq. ( IT2I) . Such deviations can be attributed to the nonlinearity, which gives rise to new 
phenomena that are not present in an ideal gas. For instance, as discussed in Sec. IIIl At for 
A(pin ^ ±7r, the merging of two independent condensates is impossible, due to the formation 
of a dark soliton. 

A pertinent question is whether the behavior of a with respect to can be explained 
solely by the movement of the traps. To answer this question, we have investigated the 
dynamics of a single atomic cloud in an adiabatically moving harmonic trap with 

V{x,t) = [x-s{t)-l]y2, (13) 

where s{t) is given by ([2]). The movement takes place from t = to t = T^, and the 
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condensate is prepared initially in the state $(x,0), which is the lowest-energy solution of 
the harmonic potential V{x, 0) = (x — /)^/2. 

For a trap moving with a constant speed v, the wavefunction of the cloud in the laboratory 



frame, is given by 



15| 

^iab(a;, t) = ^{x - ^;^)e-'('^+5^')*+^''^ (14) 



where iIj{x) is the corresponding wavefunction in the moving frame. We see therefore, that 
the main effect of the movement is to increase the chemical potential by while creating 
a phase modulation (grating) determined by vx. Both of these terms will have influence on 
the phase of the cloud at the end of the movement. 

In our model, the speed is not constant and can be defined as the rate at which the trap 
minimum changes in time, which according to Eq. ([2]) is 

v{t) = ±ls{t) = T:!^^i^(^]- (15) 



2T V T 

Given the adiabatic nature of the motion, the phase shift acquired at Tm due to the term 
can be estimated as 

A0^ovc [ v{tfdt, (16) 
^ Jo 

and using Eq. (fT5|) we obtain 

A</)movc ~ (17) 

For the range of parameters used throughout our simulations, A0move ~ 0.4, while A0move ~ 
as opposed to the linear increase of a with respect to T^. 

The analytic treatment of the phase modulation vx is far more complicated due to its tem- 
poral and spatial dependence. To investigate its role on the phase of a moving condensate, 
we have solved numerically the GPE for a condensate in the aforementioned adiabatically 
moving harmonic trap. The total phase shift that the condensate has acquired at the end of 
the movement can be read out in the way described in Sec. IIIIl and is plotted in Fig. [7] for 
various values of Tm, together with the estimation ( |T7l) . Clearly, the main role of the product 
vx is to add stepwise modifications on the power law of Eq. ( ITTl) and thus, in any case the 
linearly increasing behavior of a with respect to the merging time cannot be explained in 
the framework of moving independent condensates. 

In closing this section, we would like to note once more that the present theoretical frame- 
work cannot answer the question whether the phase of a condensate that was created by 
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FIG. 7: Phase shift experienced by a condensate trapped in a moving harmonic potential described 
by Eq. (|13p. as a function of the movement time for / = 4. (crosses) Numerical solution obtained 
by propagating the GPE (squares). An estimation given by Eq. (jl7p . 

merging two independent condensates is completely randomized or depends on the initial 
relative phase. A random new phase would be a trivial case, renders the process equivalent 
to the creation of a new condensate from a thermal cloud. Assuming that this is not the 
case, we found a relation between the initial and the final condensates. Due to our numer- 
ical treatment, in particular due to the propagation in imaginary time, results concerning 
dependence on time have to be taken with a pinch of salt. Although our aim was to perform 
the merging adiabatically, it turned out that initial phase differences alone already lead to 
nontrivial excitations in the resulting condensate. Furthermore the extraction of the phase 
iff from the numerical results turned out to be very sensitive to even small excitations of 
the condensate. Resorting to imaginary time evolution, our final result, given by Eq. (fT2|) . 
separates the final phase of a condensate created by merging two independent condensates, 
in a global phase shift depending linearly on the merging time, and the arithmetic average 
of the two intial phases. The dependence on the experimentally well controllable merging 
time places the answer to questions concerning the final phase as well as the nature of the 
dissipation processes within reach of current experiments. 
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IV. SUMMARY AND OUTLOOK 



We have investigated the adiabatic merging of two independent nearly identical Bose- 
Einstein condensates at zero temperature, within the framework of a one dimensional Gross- 
Pitaevskii equation. We have been able to answer some of the questions raised previously 
in the literature js, 7|, pertaining to the type of phase-sensitive excitations created during 
the merging, as well as the factors that determine the phase of the final condensate. 

Our simulations show that the initial phase difference between the two condensates A(^in, 
is a crucial parameter dominating an adiabatic merging process. Irrespective of the strength 
of interatomic interactions, it may prevent merging altogether, with the final density distri- 
bution exhibiting two distinct peaks. For non-interacting gases, this is due to the fact that 
the final state is a superposition of the first two eigenstates of the Hamiltonian, the excited 
one having a node at x = 0. In the interacting case, however, the separation is due to a soli- 
ton, i.e. a nonlinear phenomenon. For both cases, besides this type of excitations, we also 
have a dipole oscillation of the final condensate. All of these excitations are phase-sensitive 
in the sense that their dynamics and characteristics are determined mainly by Ayjin. Our 
estimates for the excitation energy are in good agreement with recent experimental obser- 
vations [3|, [TJ. 

Although in the one-dimensional model under consideration, phase-sensitive excitations 
have long life times, in realistic experimental setups they appear to decay quickly, increasing 
thus the temperature of the system. Removing the more energetic atoms by evaporative 
cooling, one may thus dissipate the excitation energy, without affecting the phase infor- 
mation carried by the final condensate. We simulated such a phase-information-preserving 
dissipation mechanism by propagating the Gross-Pitaevskii equation in imaginary time. In 
this way, we were able to analyze the phase of the final condensate, and show that for a 
given interaction strength it is determined mainly by the merging time, as well as the initial 
phase difference A(^in. Moreover, we derived an analytic expression for this dependence. 
The verification of this formula as well as the underlying assumptions are technically within 
reach of current experiments. 

Finally, it should be emphasized that the above results are valid over a wide range of 
interatomic interaction strengths. There are still, however, many open issues pertaining to 
the merging of two condensates with different chemical potentials, and the decay of phase- 
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sensitive excitations. These questions as well as the extension of the present one-dimensional 
model of merging to a full three-dimensional theory, are currently under investigation. 
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